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^ Abstract 

r^ The softness parameter is based on tlie relative intensity of several optical emission lines 

O emitted by the gas ionized by young massive star clusters and can be used to derive the 

^ equivalent effective temperature {T.^,) in those objects whose stellar population cannot be 

c/3 resolved. This method has several uncertainties due to the disagreement between different 

1 I synthesis model atmospheres but it is robust to study the relative variations between objects. 

Following the 2D photoionization models of the giant Hii region NGC 595 (Perez-Montero 
K^ et al. 2011) we show that the determination of T* with the rj parameter is also robust in 

T-^ different regions of a same object with large variations in the geometry of the gas and in 

^^ the dust-to-gas ratio. 

in 

p 

l>! 1 Introduction 

O 

I The study of the properties of young massive star clusters in our Local Group can be done 

;;v by means of the census of the stars that belong to the corresponding object and the analysis 

• ^H of their properties using a colour-magnitude diagram. In the case of those objects whose 

rN stellar population cannot be resolved with the facilities that are nowadays available, it is 

C^ necessary to appeal to other tecniques. A possible approach is to study the massive stars by 

inspecting their effects on the surrounding interstellar medium (ISM), including the supply of 

mechanical energy from stellar winds and supernovae explosions, the ejection of new metals 

produced in the interiors of the stars, or the excitation and ionization of the atoms in the ISM 

by the energetic photons emitted by the stars. These factors make the optical spectrum of an 

Hii region to be frequently characterized by the presence of bright prominent recombination 

lines emitted by hydrogen and helium and of collisional lines emitted by the metals. Although 

the shape of the stellar continuum and some stellar features are sometimes detectable [e.g. 

the Wolf-Rayet features), the most important source of information about the properties of 

the ionizing stellar clusters can only be the relative intensities of these lines. 
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These intensities are mainly dominated by the so-called functional parameters, includ- 
ing the metallicity (Z), the ionization parameter ([/, i.e. the ratio between ionizing photons 
and the density of particles), and the equivalent effective temperature (T*). This last is the 
unique parameter that only depends on the spectral energy distribution (SED) of the ionizing 
star or cluster. 



2 The softness parameter 

This parameter, denoted by the greek letter 77, was introduced by Vilchez &: Pagel (1988) to 
derive T^, using only the information from emission-line intensities in the optical spectrum. 
It is defined as: 

_0+/0^+ _ [0 1 1]3727 A / [01 1 1]4959, 5007 A 
'^~ S+/S^+ ~ [5//]6717,673ll/[5///]9069,9532 ^^^ ' ^' 

depending mainly on the ionic abundances of oxygen and sulphur ions. This ratio is basically 
a function of the ratio between the number of ionizing photons of O^ (LP. = 35.1 eV) and 
S+ (LP. = 23.3 eV) and has also a certain dependence on the overall metallicity [o(Z)] when 
it is derived directly from the emission-line intensities. 

In the ideal case of a blackbody it only has a linear relation with the slope of the SED, 
but for young massive star atmospheres it noticeably changes for different model atmospheres. 
In upper panels of Fig. 1 we show the SEDs of WM-Basic (Pauldrach et al., 2001 with an 
expanding spherical geometry) at left and Tlusty (Hubeny &: Lanz, 1995, with a plane-paralel 
geometry) at right for a same value of T^, (40000 K) and different values of the metallicity. 
The slope and the shape of the SEDs can be very different in the energetic range covered 
by the optical r/ parameter (Simon-Diaz &: Stasihska, 2008). An alternative to the optical 
Tj parameter was proposed for the mid-infrared spectral range with emission lines of [Nell] 
(12.8 ^m), [Nelll] (15.6 fim), [Sill] (18.7 ^m), and [SIV] (10.5 /im) (Martin-Hernandez et 
al., 2002; Morisset et al. 2004), but this involves an even more energetic range between the 
ionization potentials of S^^ (LP. = 34.8 eV) and Ne (LP. = 40.5 eV), where discrepancies 
between the different synthesis model atmospheres are larger (see Perez-Montero &: Vilchez, 
2009, for a more detailed discussion). 

In the lower panels of Fig. 1 we show the relation between the optical emission-line 
ratios of oxygen ([OII]/[OIII]) and sulphur ([SII]/[SIII]) for a sample of objects with emission- 
line like spectra compiled by Perez-Montero &: Vilchez (2009). The solid lines are quadratical 
fits to different sets of Cloudy (Ferland et al., 1998) photoionization models with different 
values of T* (from 35 kK to 50 kK). At left, these models are calculated using WM-Basic stellar 
atmospheres and, at right, Tlusty (see Perez-Montero & Vilchez to see a deeper descripion 
of these models). As can be seen, there are substantial differences between the two studied 
synthesis model atmospheres. These differences prevent an accurate determination of the 
absolute value of T^, in Hii regions but, apparently, allows a good relative characterization for 
different families of objects. In both panels, Circumnuclear Star Forming Regions (CNSFRs) 
and Hii galaxies have the higher temperatures and, in contrast, smaller regions in our Galaxy 
and the Magellanic Clouds have lower temperatures. 
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Figure 1: Upper panels: SEDs of synthesis model atmospheres from WM-Basic (left) and 
Tlusty (right) for different values of the metallicity and T* = 40 kK. The vertical solid lines 
represent the ionization potentials of S"*" and O"*". In the lower panels we show the relation 
between the ratios of optical emission line-intensities [Oil] / [OIII] and [SII] / [SIII] for different 
objects. The solid lines show the quadratical fits to sets of photoionization models with stars 
of different equivalent effective temperatures from WM-Basic (left) and Tlusty (right). 
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Figure 2: Variations of the 8 to 24 //m ratio (at left) and optical [OII]/[OIII] at right for 
different elliptical apertures around the central ionizing cluster in the giant Hii region NGC 
595. Black points represent the observed properties and the red solid line the results from 
tailored photoionization models for each of these apertures (from Perez-Montero et al., 2011) 

3 The case of NGC 595 

One ot the possible caveats when using the r] parameter to find out T* in massive star clusters 
is due to that the information comes from the ionized gas. This implies that possible variations 
in the properties of the gas other than the functional parameters, such as gas geometry or 
density, or dust structure could also affect to the r/ parameter. 

We used as a test to check to what extent these variation can affect to the determination 
of r=K the Cloudy photoionization models of the 2D structure of the giant Hii region NGC 595 
(Perez-Montero et al., 2011). The observational inputs for these models come from integral 
field spectroscopy taken with the PMAS instrument at the CAHA 3.5m telescope in the 
optical range between 3700 A - 6800 A (Relaho et al., 2010) and photometric information 
from Spitzer space observatory at 8 and 24 //m filters. In these models it was assumed a 
single ionizing source for all the Hii region according to the properties of the CMD diagram 
derived by Malumuth et al. (1996). These models reproduce the variations in the optical and 
mid-IR properties for different elliptical aperture regions with similar observational properties 
by assuming different matter-bounded geometries combined with different dust-to- gas ratios, 
in good agreement to the dust-to-gas ratio derived in the integrated region (Perez-Montero 
et al., 2011). In Fig. 2, we show the comparison between the models and the observations in 
each elliptical aperture. 



In Fig. 3 we show the relation between the two optical emission- line ratios [OII]/[OIII] 
and [SII]/[SIII] involved in the rj parameter with the fits to the photoionization models with 
Tlusty atmospheres shown in lower right panel of Fig. 1. The black circles represent the 
values derived from the tailored models for each region in NGC 595. As can be seen, with the 
exception of some models, the most part of the elliptical apertures lie in regions consistent 
with a single value of T^, across the nebula, independently on the inner variations in the 
geometry of the gas and in the dust-to-gas ratio. The value of T^ (between 35 kK and 
40 kK) is consistent to the value associated with the properties of the cluster derived by 
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Figure 3: Relation between the optical emission-line ratios [OII]/[OIII] and [SII]/[SIII]. 
Solid lines represent the same fits to grids of photoionization models of different T* for Tlusty 
atmospheres represented in Fig. 1. Black circles represent the models for NGC 595 derived 
by Perez-Montero et al. (2011) 

Malumuth using a CMD diagram, which is around 40 kK. It must be noted, however, that 
neither the slope of the grid of models for NGC 595 nor the slopes of the grids of models for 
different T* are equal to 1, which is the value of the slope for the lines with different r/ across 
the diagram, so this implies that this parameter must not be used as one-dimensional. In 
contrast, its analysis must be undertaken using a 2D approach as in this case. For instance, 
in the models of NGC 595, there is a variation in 77 in a factor 2 which is not contradictory 
with an homogeneous value of T* in the models. 



4 Conclusions 



The determination of properties of a massive young stellar cluster can be obtained by using 
information from the emission line spectrum from the ionized gas in those cases where the 
stellar population of the cluster cannot be resolved. The 77 parameter as analyzed in 2 
dimensions is a robust method to derive the equivalent effective temperature of a cluster, 
even in those cases where there are variations in the geometry of the gas or in the dust-to- 
gas ratio. The most important caveat appears as a consequence of the disagreements in the 
involved energetic ranges between different synthesis model atmospheres, but it is possible to 
use this to study the relative variations of T* in homogeneous samples of objects. This is the 
case of the variation of T,,, across spiral disks found by Perez-Montero & Vilchez (2009), which 
points to a hardening of the ionizing radiation above all in spiral galaxies of low dynamical 
mass and luminosity with late morphological type. 
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